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ABSTRACT 

We describe how to extend the excursion set peaks framework so that its pre- 
dictions of dark halo abundances and clustering can be compared directly with 
simulations. These extensions include: a halo mass definition which uses the 
TopHat filter in real space; the mean dependence of the critical density for col- 
lapse Sc on halo mass to; and the scatter around this mean value. All three of 
these are motivated by the physics of triaxial rather than spherical collapse. A 
comparison of the resulting mass function with TV-body results shows that, if one 
uses 6c{m) and its scatter as determined from simulations, then all three arc nec- 
essary ingredients for obtaining ~ 10% accuracy. E.g., assuming a constant value 
of Sc with no scatter, as motivated by the physics of spherical collapse, leads to 
many more massive halos than seen in simulations. The same model is also in 
excellent agreement with iV-body results for the linear halo bias, especially at 
the high mass end where the traditional peak-background split argument applied 
to the mass function fit is known to underpredict the measured bias by ^ 10%. 
In the excursion set language, our model is about walks centered on special po- 
sitions (peaks) in the initial conditions - we discuss what it implies for the usual 
calculation in which all walks contribute to the statistics. 
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1 INTRODUCTION 

The evolution of tlie abundance of virialized clusters is 
a powerful probe of the combined effects of the nature 
of the initial conditions (Gaussian or not?), the nature 
of gravity (general relativity or not?), and the expansion 
history of the Universe (cosmological constant or more 
complex?). This has motivated studies to understand and 
parametrize this dependence, so that datasets can more 
precisely constrain the physical models. 

There are three widespread approaches to modeling 
cluster abundances. The most accurate of these is the 
brute force numerical simulation method; one simply 
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asserts that the dark matter halos which form in these 
simulations represent clusters in the Universe. However, 
this method is computationally expensive owing to the 
nonlinear nature of gravitational structure formation. 
The other two are simpler analytic models which 
consider the statistics of the initial fluctuation field 
to model halo formation. One is based on peaks, and 
attempts to model the comoving n umber density d n/dm 
of objects of a give n mass m JSardeen et al" * 



198e 



Appel fc Jone3 . Il990l : iManriaue et all . Il998l : 1 Hanamil. 



200% . The other, known as the excursion set approach, 
seeks to count the mass fraction /(m) that is bound-up 



abundances of objects themselves (IPress & Schechtei. 


19741: Peacock & Heavensl. 19901: Bond et all. 


I199I: 


Sheth & Tormenl. 


2OO2I: iMaeeiore & Riottol. \ 


2010a: 


Musso & Shethl I2OI2I). The two are related bv the 



simple fact that mass weighting each halo must yield the 
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mass fraction in halos, so 



dm 



dn(m) 
dm 



m — p 



f{m) dm. 



(1) 



However, this subtle difference has led to rather dif- 
ferent formulations of the problem. Only recently have 
the two approaches been w ritten in the same formalism 
( Paraniape fc Shethllioij ). allowing a direct comparison 
of how they differ. Whereas all previous excursion set cal- 
culations of / were based on the stati stics of all positions 
in the initial density fluctuation field. lParaniape fc ShethI 
(|2012l ) showed that the Excursion Set Peaks (ESP) calcu- 
lation of / explicitly uses the statistics of a small subset of 
positions in the initial field (those associated with peaks). 

Working with a special subset of positions is moti- 
vated by measurements extracted from simulations show- 
ing that the excursion set prediction for the mass of the 
object in which a rando mly chosen p article ends up is 
almost always incorrect (|Whitel . 1 19961 ) , whereas the pre- 
diction for the subset of particles which are at the center 
of mass of the proto-halo is a lmost always quite accurate 
ISheth. Mo fc Tormeiil . bOOlh . Moreover, for this special 
subset of (center-of-mass) particles, the predicted mass 
depends on whether one assumes a spherical or a triaxial 
collapse - and the prediction is more accurate if one uses 
the more physically appropriate (though by no means 
perfect!) triaxial collapse model. In other words, the 
physics and the statistics are simpler, and tell a consistent 
story only if one works directly with the special subset of 
positions around which collapse occurs. More recent work 
has shown that proto-halo p ositions do indeed coincid e 
with peaks in the initial field jLudlow fc Porcianl I2OIII ') , 
and their initial velo cities are biased with resp e ct to that 
of the dark matter (E 



Ludlow fc Porcianil. bOld ) 



is expected for peaks ( Desiacaues''fc''shethl 2010l ). This 



provides additional motivation for the ESP approach. 

However, this approach raises the following techni- 
cal problem. The peaks approach uses the statistics of 
the smoothed density field, as well as its first and second 
derivatives. As a result, the rms values of these quan- 
tities play an important role. In the theory, these rms 
values are related to integrals over the initial power spec- 
trum. However, they depend on the form of the smooth- 
ing filter, and, for the TopHat smoothing filter which is 
expected to be most directly related to the physics of 
gravitational collapse, some of these integrals do not con- 
verge. For Gaussian filters there is no such problem, and 
IParaniape fc ShethI (|2012l l showed that, when expressed 
in suitably scaled (and natural) units 1/, the ESP predic- 
tion for halo abundances is in excellent agreement with 
simulations, when the measured halo abundances are also 
expressed in scaled units. Namely, if one sets 



m dn(m) 
p dlnm 



d In 



d In m 



(2) 



then the ESP I'fii') is in good agreement with the dis- 
tribution of Inv measured in simulations. However, this 
agreement is only a partial success because the relation 
between halo mass and v depends on the smoothing 



filter, and the ESP model employed a Gaussian filter 
whereas simulations use a TopHat. Therefore, when 
expressed as a function of mass, the Gaussian smoothed 
ESP prediction is not a good description of dn/dlnm in 
simulations (because the Jacobian din din m depends 
on the filtering kernel). Moreover, the Gaussian ESP 
prediction assumed that halos form through a spherical 
collapse, even though this is an incorrect description of 
the physics of halo formation (|Sheth et al.l . l200ll ). In fact, 
there is substantial sc atter e ven in t he tr ia xial collapse 



model (ISheth et al.l. I2OO1I: IShethl. I2OO9I: lOalal et all . 



I2OO8I : [Robertson et all . I2OO9I : iLudlow fc Porcianil . l201lh . 



so the latter also provides an incomplete description of 
collapse. There has been some discussio n of how to in- 
clude s tochasticity int o the p r edictions dBond fc Mver^ , 



Desiacauesl. l20oi 



cluoe s tocnasticity int o tne p r eaictions (iJaona &i 
19961: IChiueh fc Led . I2OO1I: ISheth fc Tormeil 



Corasaniti fc Achitouv 



2002 



Maggiore fc Riottol. 
201 1]: Paraniape. Lam fc Sheth 



2010b 



2OI2I : ISheth. Chan fc Scoccimarrol . I2OI2I '). "but none 
of these self-consistently combines the model for the 
stochasticity with the fact that estimates of this stochas- 
ticity, from simulations or from theory, are for the special 
subset of positions around which collapse occurs. 

Therefore, the main goals of the present work are 
to include TopHat smoothing in the ESP formalism and 
incorporate the effects of more complicated collapse mod- 
els which are seen in simulations. All this is the subject 
of Section [2l Some testable predictions which elucidate 
the relation between the ESP approach and the usual ex- 
cursion set analysis of all (rather than special) positions 
in space are discussed in Section [3] Section |4] presents 
implications for the spatial distribution of these objects 
- the ESP predictions for Lagrangian halo bias - which 
can be used to further test the model. A final section 
summarizes our results. We will show numerical results 
at redshift 2 = for two ACDM cosmologies: one with 
(fim, f^A, era, Tis, /i) — (0.3,0.7,0.9,1.0,0.7) which we de- 
note WMAPl and the other with (0.25, 0.75, 0.8, 1.0, 0.7) 
which we denote WMAP3. 



2 THE UNCONDITIONAL MASS 
FUNCTION 

2.1 Notation 

We will designate Gaussian filtered quantities with the 
subscript G. By contrast, TopHat filtered objects will 
have no subscript. The Fourier transforms of the TopHat 
and Gaussian filters are W{kR) = {3/kR)ji{kR) and 
WaikRa) = e"*" ^g/^^ respectively, with ji{y) a spheri- 
cal Bessel function. 

We will frequently need the following integrals over 
the linearly extrapolated dimensionless matter power 
spectrum A^(fc) = P{k)/2Tv^: 

TopHat-filtered variance: 



2 



dlnkA'^{k)W{kRy 



(3) 
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We will exclusively use the notation u = 5c/ (Ja where 
5c = 1.686 is the usual spherical collapse threshold. 

Gaussian-filtered spatial moments: 

ale = I A\nk/\\k)e'WG{kRGf , j > 1, (4) 



so that aid = ( (V5g)^ ) and al^ = ( (V^iSg)^ )■ 
Mixed moment: 

72 ! 



0"ln 



dlnfcA^(fc)/c^WG(A;i?G)W(fci?). (5) 



Characteristic volume: 

= V3aiG/a2G ; K = {27v)^^''rI . (6) 
These are the same as defined bv lBardeen et all (|l986l '). 
Spectral parameter: 



7 



2 

0"ln 



(-^V^^g) 



where we defined the standardised variables 



(^M> , (7) 



H = S/ao 



and 



X = —V 5g/o-2G . 



(8) 



Note that this definition of 7 is similar but not id entical 
to the corresponding one in lBardeen et al ] (|l986l l. smce 
our peak heights are defined using TopHat smoothing. 



2.2 Matching filter scales 

The technical problem which we address in this subsec- 
tion is how to ensure that the peaks in the Gaussian fil- 
tered density field Sg will have an overdensity Sth = 5c 
when smoothed with a TopHat, since essentially all mea- 
surements in simulations use TopHats only. In practice, 
we need a mapping between the scales Rth and Rg of 
the two filters. All previous work accomplishes this either 
by matching the volumes: (47r/3)i?|H = (Stt)^''^ i?G. or 
by matching the variances: 5^^^ ) — ( i5g ) ■ this sec- 
ond case, the relation between -Rth and Rg depends on 
the shape of the power spectrum. But neither of these 
conditions guarantee that a peak identified on scale Rg 
satisfies 5th = ^c- 

For this reason, we construct the mapping between 
-Rth and -Rg by finding that -Rg (at a given -Rth) for 
which { (5g|5th ) ~ 5th- This is equivalent to { Jg'Jth ) = 
^ 5th ) • This definition circumvents the technical com- 
plication arising from the variance of the Laplacian of 
5th (which involves a divergent integral). For the ACDM 
P{k) we consider in this paper, -Rg ~ 0.46-Rth with a 
mild mass-dependence which we account for, in the range 
we explore, 3 x lO"^ < m/{h-^MQ) <3 x 10^^ 

Our choice of matching between -Rth and -Rg, 
namely {5g\5) — 5, means that p{5g\5) is a Gaussian 
with mean 5 and variance ( (5g ) — s. However, {5g) ~ s 
to better than 5% and, thus, p{5g\5) ~ (5d('5g — 5) (since 
the cross-correlation between 5g and 5 is very close to 
unity). This significantly simplifies our ESP analysis for 



the following reason. In principle, the excursion set peaks 
framework also requires us to keep track of the vari- 
able X = 5' I 'J {5'^ ), where 5' = d5/ds. This scalar 
(isotropic) variable correlates with /i and x, but not with 
the ot her (anisotropic) variables used bv lBardeen et akl 
( 19861 ). Therefore, one should in principle work with the 
three-dimensional Gaussian vector (/i, x, x). However, be- 
cause p{5g\5) ~ 5o{5g — 5), it turns out that setting 
X = x yields an excellent approximation (notice that this 
relation is exact for the Gaussian filter). We have indeed 
checked that our final answers for dn/d\nm only change 
by a few percent if we switch between the approximate 
and exact treatments. In the following, we will therefore 
always set x — x and only work with the two-dimensional 
Gaussian vector {fj,,x). 



2.3 Excursion set peaks with a constant barrier 

Apart from the fact that the peak height is defined using 
TopHa t filtering, there is no form al change in the deriva- 
tion bv lParaniape fc ShethI (|2012l ) of the number density 
of excursion set peaks for a constant b arrier (and this 
deriv ation is formally the same as that in lAppel fc Jone^ . 
ll990l V In this case we get 

/ESp(i^) = VJ^ESpil^) 

= (V/K)(e-''/V^/2^) 
f°° X 

X / dx — -F'(a::)pG(a:: - 7i^; 1 - 7 ) , (9) 
Jo 

where V = m/p — iirR^/S is the Lagrangian vol- 
ume associated with the TopHat smoothing filter, F{x) 
is the peak curvatu re function from equation (A15) of 
iBardeen et all ( 1986l l and pGiy — y; E^) is a Gaussian dis- 
tri bution in y with mean y a nd variance E^. (As noted 
in IParaniape fc ShethI. |2012| . F{x) / 1 reflects the fact 
that the ESP calculation averages over a special subset 
of positions.) Since we have been careful to define 1/ and 
hence the mass using TopHat filtering, the mass function 
follows in the usual way from equation 

Figure [T] compares the ESP result for a constant 
barrier B{s) = 5c = 1.68 6 (solid red) with the fit to 
-/V-body simulations from I Tinker et al.l (|2008l 'l (dashed 
gr een). For the latter we used parameters from Table 2 
of I Tinker et al.l (|2008l ) appropriate for halos identified us- 
ing the spherical overdensity (SO) definition at 200 times 
the mean density of the universe. We used cosmological 
parameters for the WMAP l cosmology mention ed in the 
Introduction, for which the I Tinker et al.l (|2008l ) function 
gives a good fit (c.f. their Table 2). Later we will also 
show results for the WMAP3 co smology which was also 
explored by I Tinker et al.l (|2008l ) - changing cosmology 
does not affect o ur results significantly. F or comparison 
we also show the ISheth fc TormenI (|l999l ) result (short- 
dashed blue) with their {q,p) = (0 .7, 0.26) which gives a 
good fit to the MICE simulations ( Crocce et al.l . [ioiol ). 

Clearly, equation (|9| predicts far too many halos 
at all interesting scales, including the high mass end 
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0.01 : 



mixed ESP (const) 

Tinker(SOM2oo„J 

ST (MICE) 



0.001 L 







0.5 

logi(,(y=v^) 



Before explaining how this can be done, we note that, 
although equation (|10|l is a reasonable approximation to 
the mass dependence of the critical density required for 



collapse jSheth et al 



around this value Jshet 



20011'). there is substantial scatter 



Th is scatter has been 

quantified by [Robertson et al.l (j200^ 1 who approximated 
5ec bjQ 



B = <5c + 



(11) 



and found that, over the range logjQ(!/^) G (—0.3,0.9), 
the value of B was Lognormally distributed around a 
mean value of ~ (5c + 0.43a"o with variance (0.09cro). The 
physical origin of this s catter is unclear. Tr iaxial collapse 
predicts some scatter (|Sheth et al.l . l200l[ ). which is re- 
lated to the statistics of the tidal field, but the measured 
scatter is considerably larger. In any case, the physical 
origin of this scatter is unimportant for what follows - 
though we believe that this is an issue which deserves 
further investigation. 



Figure 1. Excursion set peaks mass function for the mixed fil- 
tering approach discussed in the text, using a constant barrier 
(solid red). This is compared with the fit to A^-body simu- 
lations given by [Tinker et al] l|2008l '). appropriate for spheri- 
cal overdensity halos that enclose 200 times the mean density 
of the universe (dashed gre en). For comparison, we also show 
the lSheth fc TormenI lll999l ) mass function (short-dashed blue) 
with their {q,p) = (0.7,0.26), which g ives a good fit to the 
MICE simulations llCrocce et al.l |2010|) . The ESP+constant 
barrier mass function overpredicts the abundance at all scales 
of interest. All curves used the WMAPl cosmology - we have 
checked that changing the cosmology does not significantly 
alter the comparison. 



(i/'^ = 10 corresponds to m = 1.3 x IO^^H'^Mq for this 
cosmology) . 



2.4 Beyond the spherical-collapse constant 
barrier 

It is tempting to attribute the mismatch shown in Fig- 
ure[T]to the effects of non-spherical collapse. While these 
effects are certainly important at smaller mass scales, 
it is not obvious that they also matter at the largest 
scales, where the spherical collapse approximation should 
be quite accurate. To see that they are in fact relevant, 
note that the triaxial collapse model predicts that the 
critical density for collapse scales approximately as 

<5ec ~<5c[l + 0.4(ao/5c)'-^] (10) 

( Sheth et al.l . 1200 ll ). Simulations indeed show that this 
expression for S^c provides a reasonable approximation to 
the mass dependence of the critical density required for 
collapse ( Sheth et al.l . lioOll ). However, this implies that 
when V = Sc/ao is as large as ~ 4, then departures from 
the spherical Sc are already of the order of ten percent. 
Therefore, we must include deviations from the spherical 
coll apse barrier in orde r to have a fair comparison with 
the [Tinker et al.l ( 20081 ) simulations. 



2.5 Excursion set peaks with a square-root 
barrier and scatter 

To include the mean trend a nd sc atter in barrier 
hei ghts shown by ISheth et al.l (|200ll ) and quantified 
by [Robertson et al. I (I2OO9I). we will adopt the following 
model, first discussed in Paranjape et al.l ( 2012f ). We con- 
sider a family of deterministic "square-root" barriers B 
given by equation (|lip . with /3 a stochastic variable. 
Therefore, if /esp(!^|/3) denotes the ESP solution for a 
given value of /3, then our full solution will be 



/esp(/^)= / d/3p(/3)/ESp(/^|/3). 



(12) 



Below we use simulations to motivate our choice of p(/3), 
and provide an explicit expression for fBSpi^lP)- Still, it 
is easy to understand qualitatively what happens. 

For excursion set peaks, the square root barrier (jlip 
with /3 > drama tically lowers the halo cou nts at large 
masses (Figure 3 in lParaniape fc Shethl[2012l ). While this 
goes in the right direction, introducing scatter in the 
value of P will counteract this decrease. This is because a 
scatter in /3 induces a scatter in the predicted mass of the 
halo. The steepness of the mass function implies that the 
scatter preferentially moves small mass objects to larger 
masses (the same effect that leads to Eddington bias), 
which in turn increases the abundance at large masses. 

The importance of the scatter depends on the dis- 
tribution of p. Our choice for p(/3) is particularly simple. 
As discussed in the previous section, the distribution of 



^ We u se the results of [ Robertson et al. I Hooi) rather than 
those of lPalal et al.l l|2008l ) because the latter only show a plot 
whereas Robertson et al. provide a fit to a probability distribu- 
tion. Both groups find a similar mean and r.m.s. for the barrier 
heights, however, which increase approximately linearly with 

TO- 
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Figure 2. Ratio of the ESP mass function for a square-root 
barrier with scatte r in the barrier slope (equation I12|l to the 
[Tinker et all bOOSl) fit, for WMAPl (solid red) and WMAP3 
(dotted red) cosmologies. In each case we set Sc = 1.686 and 
p(/3) to be Lognormal with (/3) = 0.43 and v'Var(/3) = 0.3, 
for the reasons discussed in the text. For comparison we also 
show the ratio of the MICE fit with the Tinker et al. fit for 
the WMAPl cosmology (short-dashed blue). Horizontal dot- 
ted lines mark 10% deviations from the Tinker et al. fit, which 
is the typical scatter in their measurements in this mass range. 



B seen in simulations for halos of a fixed mass is Lognor- 
mal, so we choose the distribution of p to be Lognormal 
with {P} = 0.43 and Var(/3) = 0.09. 

For square-root barri ers with fixed /3 (equation 19 of 
IParaniape fc Shethl . [ioil . for the barrier in equation 1111 
above) , 



/espHP) = (^/K.)(e^(■^+«'/V^^) 
' da; !^^l^F{x) 



/37 



71/ 



X pg{x - /J7 - 71/; 1 - 7 ) . (13) 

We have assumed /3 > 0, which is true for the Lognormal 
model we will consider. If /? is allowed to take negative 
values as well, then the lower limit of the integral over 
X must be replaced by Xmin = max{0,/37}, so as to en- 
force the peak constraint. While there is no closed form 
expression for /esp(j^), the integrals involved can be eas- 
ily computed numerically. Note that, because we use the 
distribution p(/3) that is seen in simulations, our model 
has no free parameters. 

Figure [2] shows the results for the WMAPl (solid 
red) and WMAP3 (dotted red) cosmologies. Note that 
the value of erg is ~ 10% smaller in the WMAP3 cos- 
mology. We plot the ratio of the ESP prediction in equa- 
tion (fT2)) and the [Tinker et all (|2008l ) fit. For comparison 
we also indicate the ratio of the MICE fit with the Tinker 



et al. fit (for the WMAPl cosmology) as the short-dashed 
blue curve. The horizontal dotted lines mark 10% devia- 
tions from the Tinker et al. fit, which is the typical scatter 
in their measurement s in this mass range (c.f. figure 6a 
of [Tinker et"al] . l2008l ). This shows that our model, which 
has no free parameters, works rather well. 



3 IMPLICATIONS FOR AVERAGES OVER 
ALL, RATHER THAN A SPECIAL 
SUBSET OF POSITIONS 

The ESP model of the previous section is the first to 
self-consistently integrate the complex physics of non- 
spherical collapse with the fact that the physics is almost 
certainly simplest around the center-of-mass of the proto- 
halo. Although the agreement with the halo counts rep- 
resents a significant and nontrivial success, and the next 
section shows that the predicted spatial distribution is 
also accurate, it makes a number of other testable pre- 
dictions which are related to the logic of the excursion 
set approach. 

For example, one can construct the center-of-mass 
random walk associated with each halo and explicitly 
check, object by object, whether or not the correspond- 
ing barrier was crossed at a larger scale. This is, in 
fact, precise l y the kind of test which was performed by 
ISheth et al.l (|200ll ). Revisiting this test is left to future 
work. 



3.1 Distribution of S for all particles in 
proto-halos of a given mass 

We remarked in the Introduction that our ESP model 
explicitly averages only over special positions in the uni- 
verse (see text following equation[9]for exactly where this 
happens). If we treat these positions as the fundamen- 
tal quantities (for which the physics of collapse is most 
easily applied), then by accounting for their initial den- 
sity profiles, it is straightforward to estimate what would 
have happened had we averaged over all positions. Be- 
cause our model has been calibrated for center of mass 
positions solely, predicting the right distribution of 5 val- 
ues within each proto-halo provides another test of our 
approach. We sketch the idea below. 

We begin with noting that the distribution of S at 
distance r from a peak is a Gaussian with mean approxi- 
mately given hy a {u — ■y {x\h'))/{l — 7^) C('")/C(0) and 
variance approximate ly \l — (£(r)/£(0) )^/(l — 7^)] 
(equations 7.8-7.10 in iBardeen et all . [1989 ). Note that 
the mean value asymptotes to crv ^{r) / (,{0) in the limit 
iy» I. Since ^(r) < ^(0), the term C(r)/^(0) = H is the 
ratio by which the density at r drops from its value at 
the center. 

In our model, at any given mass scale there is a distri- 
bution of 5 values for the special positions around which 
collapse occurs; call these values (5pk. This distribution 
has a variance, say, cr'^E|. And, for each 5pk, there is a 
Gaussian distribution of S for any position at a distance 
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r < R from the center (the one just described above). 
Since aU these have the same variance, the result of aver- 
aging over the r-shells yields a distribution with variance 
times 1 - hV(1 - 7^) + 2^ ^j- 
The shape of this distribution should still be approx- 
imately Gaussian if ^ 1 — H^/(l — 7^). This should 
indeed be the case; we know that 7 ~ 1/2 and E/3 ~ 1/3, 
and we can reasonably assume 1/2 < H < 1, with H not 
being very close to unity (since we average over r-shells 
that are within one smoothing radius of the halo center 
of mass, and the shells near r — contribute little to the 
total volume). Consequently, the expected distribution of 
S, measured by looking at all positions within protohalos 
of mass m{a^) (i.e. not only those corresponding to the 
proto-halo centers of mass) is a Gaussian whose mean 
value is a factor 1/2 < H < 1 smaller than that associ- 
ated with the physics of collapse, and whose variance is 
~ (T^ (1 — 1.25 H^). For example, if H = 3/4 then the vari- 
ance should be ~ 0.3a^ . This dependence on a, as well 
as the fact that H appears both in the expression for the 
mean S as well as its variance, are both testable, as is the 
prediction that the shape should be well approximated 
by a Gaussian, even though the distribution of 5pk is not. 



3.2 Distribution of predicted mass for all 

particles in proto-halos of a given mass 

The lower mean value implies that the mass which the 
excursion set approach predicts for the particles which 
are not at the center of mass should be systematically 
smaller than the mass of the halo in which they end up 
- in qualitative a greement with one of the findings of 
ISheth et all (|200ll . see their Figures 2-4). 

To see that this is indeed the case, we can use our 
ESP approach (walks centered on peaks in the initial 
field) to model what the usual excursion set estimate 
(walks centered on all positions in the initial field) re- 
ally represents. Recall that, for the statistics of all po- 
sitions, we know that if Sc is independent of m then 
the excursion set prediction for the mass fraction in ob- 



where 



jects with mass m(s) is s/(s) = u exp{—i/ /2)/\'2tv with 
(Strictly sp eaking, to compare most directly 



with lSheth et all (|200ll ') we should account for the distri- 
bution in 5c values at each s; but the constant determinis- 
tic Sc discussion below captures the gist of the argument. 
Also, our expressions below are appropriate for a sharp-fc 
filter, but the logic, and so the qualitative trends, do not 
depend on this choice.) 

In our ESP model, this fraction is given by taking 
each peak of mass M, then accounting for the distribu- 
tion of excursion set predictions for s(m) > S{M) which 
come from all the other mass elements which make up 
M (rather than just the one mass element at the center 
of mass), and integrating the distribution of S over the 
range < S" < s. I.e., 



5'/esp('S') 



M dnESp(M) 
"i dluM 



dluAf 



dln5 



(15) 



(from equation [2I1 and ,f{s\S) is the excursion set pre- 
diction for the mass fraction of Af which is incorrectly 
predicted to instead have mass m — M (s) instead of 
M{S). Because we know sf{s) as well as nBSp('S') we can 
solve (by numerical back substitution) for f{s\S). How- 
ever, we can get a simple, intuitive estimate of the result 
as follows. 



For constant 



Figures 



and 



IParaniape fc ShethI (|2012h shows that 5'/esp(5') 
is rather well-approximated simply by vf{j^) = 
V exp(— j^^/2)/-\/27r with = i^Saf' / S for some 

H < 1. (Strictly speaking, they used the Sheth- Tormen 
functional form for vf{v); however, at high masses, 
this form reduces to the simpler one given here with 
H = VO.7.) Since this is the same functional form for S 
that appears on the left hand side of equation (|14p for s, 
it is easy to check that sf{s\S) should also be the same 
function of v, but with = {5c — H5c)^/(s — S). 

Thus, this calculation shows explicitly that the av- 
erage over all walks will yield / which has a similar func- 
tional form to the one for the specially chosen subset, 
but the two will have different values of 5c- Moreover, 
the distribution of predicted masses will be very different 
from a delta function centered on M. Rather, it will look 
just like the conditional excursion set prediction for the 
fraction of all walks which first cross 5c on scale s given 
that they are known to pass through S5c < 5c on scale 
S before first crossin g 5c- This latter po int quantifies the 



AS 



ShsHS) sfis\S). 



(14) 



betore nrst crossm g dc. inis latter po int quantines tne 
explanation given bv lSheth et al. ([gOOlh when addressing 
the concerns raised bvlWhitel lll996l)~ It also ex pands on a 
point first made bv lParaniape fc ShethI (|2012l l: the ESP 
framework naturally accounts for the fact that one must 
rescale v — y/^v with q ~ 0.7 if one wants to use the 
usual excursion set expressions (based on the statistics of 
all positions) to model halo counts - but this rescaling 
is not necessary if one uses the statistics of the special 
positions for which the physics is simplest. Of course, be- 
cause the physics of collapse refers to collapse around the 
halo centers of masses, and not around random positions 
in space, this discussion highlights the dangers of inter- 
preting the barrier height in the random walk calculation 
with that for the physics of collapse around special posi- 
tions in space. 



4 HALO BIAS 

In this section we derive expressions for the predicted 
clustering of halos in our model. 

The Lagrangian bias associated with excursion 
set peaks follows from writ ing the con di tional mass 
function as discu s sed b y IMusso et all (I2OI2I ) and 
IParaniape fc ShethI (|2012l '). Introducing scatter in the 
barrier parameter /3 is straightforward, as we discuss 
below. The ESP conditional multiplicity function for a 
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square-root barrier at fixed /3 and a fixed density i5o at 
scale 5*0 is given by 

hspW,5o) = {V/V,)p{B/V^s\5o) 
1 



(16) 



X — / dx(x ~ l3j)F(x) 
X p{x\B/^, So) 



wliere p(/x|5o) and p{x\^.,5o) are 1-dimensional condi- 
tional Gaussian distributions tfiat can be written us- 
ing the covariance matrix of the trivariate Gaussian in 

The arguments of lMusso et alj ( 20121 ) show that the 
bias coefficients are the Taylor coefficients of the 

expansion of equation (|16p in terms of 5q = So/S^, upon 
imposing the condition c = 0, where c is the matrix that 
one must subtract from the unconditional covariance ma- 
trix of the variables (/^, x) to obtain the covariance matrix 
of the conditional Gaussian x-|5o)- The resulting con- 
ditional multiplicity is 

/esp(/^|/^,5o,c = 0) 

= (l^/K)(e-(-+''-'"''(^x/^'''''''/V%^) 
1 f°° 

X — / da; (a; — l3'y)F{x) 

X pg{x ~ Pj - -^ui;! - j^) , (17) 

where we have defined ui = u{l — So{Sx / So)(l — ex)), 
and 

Sx= j dlnkA'^{k)W{kR)W{kRo), 

ex = 2dln5'x/dlns (18) 

are cross-correlations between the small and the large 
scale. 

At fixed /3, the bias parameters can be written as 



(Sx/So) 
(Sx/5'o)2 



= ^^l{u,^3) + {l-ex)Xl(u,^3), (19) 
= fi2{u,l3) + 2(1 - ex)/ii(i',/?)Ai(/^,/?) 



+ (1 



'A2(/^,/3), 



(20) 



where, using F = 7/ \/l — 7^ 



lin(u,l3)^u''H4iy + l3), (21) 
= {-VuYiHr. {y-PT-Ti.)\u,p)^ , (22) 



where Hn{x) 



e''''^^{-d/dx)"e- 



V2 



are the "prob- 



abilist's" Hermite polynomials, and for some function 
h{y, v, /?), we have 

(%,i.,/3)ii.,/?)^ 

/°; dy (y - pr)F{y^/r)pG{y - /3r - Tu; l)hiy, v, /3) 



j;^ dy {y - pV)F{y^/V)pG{y ~ PT - Tu- 1) 



15 
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Figure 3. Scale-independent (peak-background split) linear 
Lagrangian halo bias <5c6io in the ESP framework for a square- 
root barrier with Lognormal scatter (solid red). This curve is 
in excellent agreem ent with the fit to A'^-body simulations from 
[Tinker et"al] I2OI0I '). the simulations being the same as those 
for which Figure [2] showed a mass function comparison (see 
text for a discussion) . Th e dotted black curve is the spherical 
collapse prediction u'^ — 1 \Mo fc Whitd . Il996l) , and the short- 
dashed blue curve is the peak-background split prediction as- 
sociated with the mass function fit to the MICE simulations. 



Marginalising over j3 gives the bias parameters as 

'Sx 



5.b,{v) = ( ^ 



{^ii\v) + {l-e,){M\v) 



(24) 



(/i2li^) + 2(l-ex)(AiiAi|i/) 



+ (l-ex)'(A2|i.) , (25) 

where, for some function g{v, /3), we have 

, , , Jd/3p(/3)/ESp(^|/?)g(^,/3) 
\9\vi^ CAR m\f rim ■ 

The presence of ex in these expressions for real-space 
bias is an indication of fc-dependent Fourier -space bias, 
as d iscussed at length bv lMusso erahl (|2012l ). 

IMusso et all (|2012l ) also showed that the usual peak- 
background split expressions for bias (which they de- 
noted as b„o) correspond to taking the large scale limit 
fJo — >■ oo, or equivalently, setting ex — > and Sx/Sq — 1 
in the coefficients b„. Figure [3] shows our ESP prediction 
for the linear peak-background split bias (5c6io (i.e., set- 
ting ex = 0, Sx/So = 1 in equation [24]) as the solid red 
curve. The dotte d black curve is th e spherical collapse 
prediction - 1 (|Mo fc Whit^ - llOOd). The da s hed g reen 



curve shows the fit presented bv [Tinker et al. (l201ol) fo r 



(23) 



the same simulations analysed bv [Tinker et al. ( 20081 ) . 
We have checked that in this case the results for the 
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Figure 4. Same as previous Figure, but now for the second 
order bias coefficient, fn this case, there are no measurements 
from simulations with which to compare. 



WMAPl and WMAP3 cosmologies are nearly identical, 
so we only show the results for WMAP3. The ESP curve 
matches the A'^-body fit very well at large masses, which is 
where the ESP fra mework is expected to work well. Note 
that the curve for iTinker et all (|20ld ) is not simply the 
logarithmic derivat ive of the mass function calibrated by 
[Tinker et al.l (|2008l ) ; that derivative would underes t imate 
the large sca l e bias by up to 10% (|Manera et al.1 . 12OI0I : 
[Tinker et all I2OI0I ). (For comparison, the short-dashed 
blue curve shows the peak-background split prediction 
for the MICE mass function.) 

Figure [4] shows the corresponding result for the 
quadratic peak-background split coefficient ^^^20 for ESP 
(solid red) and the spher ical collapse result v^iy^ — 3) 
( Mo. Jing fc Whitel.ll997l. dotted black). The dashed blue 



Mo et all j 19971 ) expression after rescal- 



curve shows the 
ing V — >■ V0.7i/. In this case, there are no measurements in 
simulations with which to compare - this is the subject of 
ongoing work - but we have included the Figure to show 
that &20 behaves qualitatively like 610, transitioning from 
the usual (sharp-fc excursion set) expression rescaled by 
V — > VO.7 1' at ~ 1 to lying above this rescaled value at 
larger v. 

These pl ots illustrate another point made by 
IParaniape fcShcth ((201^ 1: the ESP framework naturally 
accounts for the fact that the measured large scale bias 
deviates from that predicted by a naive application of 
the peak-background split to a mass function that was 
fit using the rescaled variable v — > yfqv with q ~ 0.7. 



4.1 Reconstructing the peak-background split 
parameters 

In this subsection we revisit the idea of reconstructing the 
peak-background split parameters 6„o from suitable real- 
space measurements of the parameters 6„ (|Musso et al.l . 
I2OI2I ). In particular, we wish to address some subtleties 
that arise due to the stochasticity of the barrier height. 

The theoretical quantities plotted in Figures [3] and [4] 
required us to set tx — >■ and Sx/So — ^ 1 by hand. In 
Fourier space, this would correspond in practice to taking 
the limit k ~ I/Rq 0, which is how lTinker et al.1 (|2010|) 
e.g. d erived their fit for the linear bias. Musso et al.l 



(I2OI2I ). on the other hand, pointed out that the theo- 
retical predictions for S^b„ at finite Rq are equivalent to 
an average over the Hermite-transformed initial density 
field Hn{5o/^fSo) at the centers of halos. To be precise, 
if there are A'^ halos in a given mass bin, one estimates 
the bias parameters at this mass as 



, (msd) 



= s-„ 



n/2 



1 ^ 



N 



(27) 



where 5oi is the Lagrangian density contrast smoothed 
on scale Ro and centered on the i"^ halo in the bin, and 
5*0 is the linearly extrapolated variance at Rq. 

If we use a deterministic barrier in our theoreti- 
cal model for the bias, then these measurements can 
be explicitly converted into esti mates of the peak- 
background split parame ters fono (|Musso et al.l . I2OI2I : 
IParaniape fc ShethL l2012h . To see this, note that for 
the square-root barrier with fixed /3, we can use equa- 
tions (|19|) and (|20|) to write the peak-background split 
parameters in terms of 61(1^, /3) and 62 (j', P) as 



5cWQ{v,li) = /^i + Ai 
1 



(1-^x) 



(5x/So) 



(28) 



Slb2o{u,P) = fj.2+ 2^iiAi + A2 



(l"^x)2 



(Sx/S'o)2 



-{2 - ex)^2 



(29) 



which red uce to the expres sions in equations (44) 
and (47) of lMusso etahl (|2012l ) when /3 = 0. 

If /3 is stochastic (as we have seen it must be), then 
the measurements on the r.h.s. of equation (|27|l corre- 
spond to the marginalised bias coefficients - bi^^"^^ o 
(bnif, l3)\i/} - since we did not keep track of the value of 
P for each halo. For the linear coefficient 610 this does not 
pose any problem, and we can estimate this coefficient us- 
ing {nilu), which is calculable, and the measured linear 
bias. 



Scb 



(msd) 



(l-^x) 



(msd) 



(5x/So) 



(30) 
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For 620, however, we must deal with a term involv- 
ing {bi{i', /3)iJ,i{i', j3)\i') . In practice, measuring this term 
would involve knowing the value of jS for each halo. 
While this is do-able in an A'^-body simulation, it 
makes the reconstruction technique very cumbersome. 
In order to keep the algorithm simple, we will assume 
that this term can be split as {bi{v, I3)fii{v, I3)\p) = 
{bi{i', I3)\i/) { ij,i{v, jj)]!/) = b'^^^'^'' (nilv). This assump- 
tion can be explicitly tested using A''-body simulations. 
With this caveat, the reconstructed quadratic coefRcient 
is 



r2 7 (msd) 
Or Don 



2e>, 



2, (msd) 



Stb^ 



(5x/So)2 



4&i 



(msd) 



7^ (Mik) -{fJ-iW) 
/oo) 



ex(2 - ex) {/i2|!^) 



(31) 

The remaining averages over /3 simplify to some extent, 
and we have 

{v^ -l-2u{P\u) + {p^\v)) , (32) 



with 



Jd/3p(/3)/ESp(t^|/3)/3^' 
/d/3p(/?)/ESp(i^i/3) 



(33) 



5 DISCUSSION 

We have presented a derivation of equation (|12[) . which 
represents the first analytic model of halo abundances 
which accounts self-consistently for the fact that the 
physics of collapse is best described around the halo cen- 
ter of mass, is not spherical, and can vary substantially 
from one place in the universe to another. The compari- 
son between equation (|12|l and halo counts in simulations 
(Figure [2| is very encouraging. (Figure [1] shows that ac- 
counting for the physics of non-spherical collapse is nec- 
essary.) 

We also derived expressions for the associated large 
scale bias factors &io and 620 (these follow from setting 
Ex — >■ and Sx/So — s> 1 in equations 1241 and 125^ : the 
former is in very good agreement with measurements 
from simulations. We are in the process of determin- 
ing if 620 is similarly accurate. Additionally, we showed 
how real-space measurements of bi and b2 can be con- 
verted to practical estimates of 610 a nd &20 (equat i ons | 3pi 
and 131 p. which extends the res ults of lMusso et al.l ( 20121 ) 
and IParaniape fc ShethI (|2012l ') to the case of stochastic 
barrier heights. Testing the accuracy of this reconstruc- 
tion in A'^-body simulations is also work in progress. 

Our ESP analysis highlighted the benefits of devel- 
oping a model for the abundance of positions in space 



around which the physics of collapse is simplest. In Sec- 
tionO we showed that it can be used to understand what 
would have happened if one had instead studied all initial 
positions. This provides the basis for a number of other 
tests of the more formal aspects of our approach - as well 
as elucidating the relation between our ESP calculation 
and the usual excursion set approach. 

Our model has no free parameters, although it must 
take as input one physical quantity from simulations. 
This quantity describes how the typical density required 
for collapse depends on proto-halo mass, and how much 
this density varies from one proto-halo to the next. While 
the triaxial collapse model correctly predicts that this 
mean value and the scatter around it should increase as 
halo mass decreases, the predicted scatter is smaller than 
observed. We therefore hope that our work will motivate 
studies of the physical origin of p(/3). One obvious di- 
rection here is that peaks have non-trivial deformation 
and inertia tensors. Non-spherical initial shapes mean the 
critical density for collapse will now depend not just on 
the tidal field but on the shape and the misalignment 
angle between the principal axes of the two tensors. In 
principle, our ESP model contains a prescription for the 
distribution of shapes (and misalignments) for which one 
should run the ellipsoidal collapse model. 

Futhermore, we do not expect our model to apply 
at masses for which 5c{m)/ao{m) < 1, as there is no 
physically compelling reason why a peak of (normalized) 
height 1 should dominate its surroundings. While this 
mass regime may be less interesting for studies which use 
the abundance of rich clusters to constrain cosmological 
models, understanding it is useful for understanding the 
formation and evolution of galax;ies. So we hope our ap- 
proach of identifying just what it is that is special about 
the positions around which gravitational collapse occurs, 
and then using only these positions to make inferences 
about halo assembly histories, will eventually provide in- 
sight into the low mass regime as well. 
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